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Abstract 

The problem of path generation for the spherical 4R mechanism is solved using the Differential 
Evolution algorithm (DE). Formulas for the spherical geodesies are employed in order to obtain the 
parametric equation for the generated trajectory. Direct optimization of the objective function gives 
the solution to the path generation task without prescribed timing. Therefore, there is no need to 
separate this task into two stages to make the optimization. Moreover, the order defect problem can 
be solved without difficulty by means of manipulations of the individuals in the DE algorithm. Two 
examples of optimum synthesis showing the simplicity and effectiveness of this approach are included. 

1 Introduction 

Spherical mechanisms have a large number of applications (see for example, [1-10]), being the spherical 
4R mechanism one of the most studied. Nevertheless, the vast majority of these studies deal with function 
generation, and there are relatively few studies on the prescribed path generation task for optimal synthesis 
of spherical mechanisms [11]. The optimal synthesis of spherical mechanisms for the path generation task 
without prescribed timing, has been addressed in [12]. The synthesis procedure presented there consists 
of two layers of optimization, which are handled separately. In an inner layer, the objective function is 
minimized over the input angles with fixed linkage dimensions for a given iteration of the outer layer. The 
output of this procedure is the set of input angles defining the linkage configurations under which the 
coupler point attains the closest position from a given point. On the other hand, the outer layer minimizes 
the linkage error over the design vector. In this process, iterations do not stop until some conditions are 
met. Recently, two studies have been published regarding the prescribed timing case [11,13]. In these 
studies, the Fourier coefficients are used to describe the coupler curve, then the path generation task is 
done by using the atlas method. 

In this study we address the problem by optimizing the structural error which is defined as the sum 
of the square of the distances between the desired and the generated points. The generated points are 
obtained by finding the spherical geodesic between two points. The optimization is done by using the 
evolutionary algorithm known as Differential Evolution (DE) , that allows us to modify the main operators 
and thus provide a simple way to deal with the order defect problem. 

When dealing with the task of path generation without prescribed timing, the order defect problem 
has to be solved. The methodology to obtain the parametric equation for the generated trajectory that 
is proposed here, allows the designer to deal with this problem without difficulties. The optimization 
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process can be done by minimizing a single objective function without the need of making a sequential 
optimization. 

The paper shows in detail how to construct the objective function (fob) for optimal synthesis of 
spherical mechanisms in both cases: the prescribed timing and the not prescribed timing path generation. 
We also work out two examples of optimum synthesis by using the proposed methodology. 



2 Mathematical preliminaries 

In this section, relevant concepts and useful formulas involved in the synthesis of spherical mechanisms 
are presented. In particular, a rotation matrix is used to construct the parametric equation for spherical 
geodesies which makes it possible to construct the parametric equation of the generated trajectory. 



2.1 The rotation matrix 

In the study of synthesis of a spherical four bar mechanism one kind of transformation is very important: 
the one that keeps the the length of a vector invariant. It is well known that rotations of the coordinate 
system, or rotations of a vector, keep the vector length invariant. Both transformations are carried out 
by orthogonal matrices (A T A = 1). If the coordinate system is rotated, the rotation is known as passive; 
if the rotation is performed on the vector, the term active rotation is used [14]. In the passive case, the 
vector remains stationary and rotations are made on the coordinate axes. In the active case, the vector 
is rotated while the coordinate system remains fixed. Relationship between them is given in Eq. (2.5). 

In order to find the rotation matrix, the generators of the SO (3) group (the set of 3 x 3 orthogonal 
matrices with determinant equal to 1) are required. 

The generators are obtained from the formula: 
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Let us consider a passive rotation of an angle 9 around a unit vector n = [n x n y n z ] . The rotation 
matrix is then written in terms of the generators as: 

R passivc (#,h) =exp (*0$^T fc n fc ) . (2.4) 

V k=x ) 

On the other hand the active rotation would be: 

R(0,n) = R passivc (-0,n). (2.5) 

It is worthwhile to mention that there are other forms to describe the rotation of a vector, an interesting 
discussion can be found in [15]. 

2.2 Spherical geodesies 

It is assumed that with a point (x, y, z) on the space there is a vector associated to it, namely, the vector 
from the origin to the point (x, y, z). As it is customary, we will use the same symbol to represent both, 
it will be clear from the context if we are talking about the point or the vector. 

In order to construct the geodesies on the spherical surface, let us consider two points hi and h 2 on 
the surface of a sphere of unit radius. The parametric equation for the spherical geodesic r gl2 from hi to 
h 2 is constructed by rotating the vector hi an angle cos -1 (hi • h 2 ) around of the normalized cross product 
vector hi x h 2 denoted by n?ji2- 

r s i 2 (0;hi,h 2 ) = R(0,n W2 )h 1 , (2.6) 
where 9 is the parameter of the trajectory varying from 9 = up to 9 = cos _1 (hi • h 2 ). 



3 Kinematic modeling 

Let Xi, x 2 , x 3 and X4 be four arbitrary points on the spherical surface. We will consider the input link 
as the geodesic connecting the points xi and x 2 , the coupler link as the geodesic connecting the points 
x 2 and X3, the output link as the geodesic connecting the points X3 and X4, finally the fixed link will be 
the geodesic connecting x 4 and xi. Assuming that we are in a unit sphere, the lengths of the links will 
be the angles between the vectors defining the links. We will call «i, a 2 , a 3 , a 4 the lengths of the input 
link, the coupler link, the output link and the fixed link respectively. Fig. 1 shows the aforementioned 
variables. 

Now, let us consider that the input link is rotated from the assembly position by an angle 9. The 
parametric trajectory for the extreme of such a link, r 2 (9; xi, x 2 ), will be obtained by rotating the vector 
x 2 around the vector xi an angle 9, 

r 2 (0;xi,x 2 ) = R(0,xi)x 2 . (3- 1 ) 

Such a rotation will generate in the output link a rotation given by the angle <j>(0). The parametric 
trajectory for the extreme of the output link is written as 

r 3 (6»;xi,x 2 ,x 3 ,x 4 ) = r 3 (0(0);x), with x = (xi, x 2 , x 3 , x 4 ). 

It is given by 

r 3 (0(0);x)=R(</>(0),X4)x 3 . (3.2) 

The dependence of equation (3.2) on x 3 and x 4 is clear, the dependence on xi and x 2 comes from (f)(9) 
(in fact, to be precise we should write </>(#; x), although it is not written to avoid unnecessary notation). 

Notice that the parameter 9 is not the angle between the geodesic associated to the fixed link and 
the geodesic associated to the input link but the rotation angle of the x 2 vector around the xi vector. 
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Figure 1: Variables involved in the Spherical four-bar mechanism. 



Similarly <j>(9) is not the angle between the geodesies associated to the fixed link and the output link but 
the rotation angle of the x 3 vector around the x 4 vector. 

The angle <j>(0) is obtained by requiring the coupler link to have a constant length. In other words, 

r 2 (0;xi,x 2 ) • r 3 (0(#);x) = x 2 • x 3 = constant. (3.3) 

Solutions for <j>(9) from Eq. (3.3) can be obtained either numerically or analytically. Indeed, the analytical 
solution can be found in [16], which in our notation is given by: 

4> = $ - 2 tan" 1 ^ ±V ^+| 2 — ^ where (3.4) 
A = sin a\ sin a 3 sin(# + O ) 

B = cos ot\ sin a 3 sin a 4 — sin a\ sin a 3 cos a 4 cos(6* + <do) 
C = sinai cos a 3 sina4 cos(6> + 0o) + 
cos at\ cos a 3 cos a 4 — cos a 2 , 

with Go being the initial angle between the geodesies associated to the fixed link and the input link. <I>o 
is the angle between the geodesies associated to the fixed link and the output link. 

In order to obtain the parametric equation for the generated trajectory r gen (as it is shown in Fig. 
2), it is necessary to find the position vector v cp {9, v\ x) of an arbitrary point on the coupler link (or on 
the complete arc that contains the coupler link, since < v < 2n). This can be done by rotating the 
vector r 2 (6*;xi,x 2 ) by an angle v around the unit vector n 23 orthogonal to r 2 (0;xi,x 2 ) and r 3 (^>(^);x). 
Eq. (2.6) then yields: 

r cp (6>,^;x) = R(z/,n 23 )r 2 (6>;xi,x 2 ), (3.5) 

where 

- = r 2 (fl;xi,x 2 ) x r 3 (^(6>);x) 

1123 ||r 2 (0;x 1 ,x 2 ) xr 3 (<X0);x)|r { ' 

Using Eq. (3.5), the parametric equation for the vector r gcn is obtained by rotating the vector r cp (8,/3 + 
7;x) an angle n/2 around the vector r cp (6, (3; x): 

r gcn (9, (3, 7 ; x) = R (tt/2, r cp (9, /3; x)) r cp {8, + 7 ; x). (3.7) 
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Figure 2: Vectors involved to obtain r gon . The dependence on 9 and x is not written. 



Since the point (k = 1, 2, 3, 4) is on the unit sphere, two parameters are required to determine its 
location. In this work the azimuthal <p (latitude) and polar rj (colatitude) angles are used to generate the 
point, hence x fe = x k (ip 7 i]). Explicitly, 

Xk(<fk,rik) = (cosLp k s'mf] k ,sm(p k smT] k7 cosf] k ). (3.8) 

In the case of the prescribed timing situation, Eq. (3.7) depends on eleven parameters, or on ten param- 
eters if all input angles are given 1 . Those parameters will be the adjustment parameters (or independent 
variables) for the optimum synthesis problem. 

4 Path generation 

The task of function generation for the spherical mechanism has been widely studied. However there 
are relatively few 2 articles addressing the task of path generation without prescribed timing. Here, we 
present how to construct the objective function for this case. The optimization of such a function is 
performed directly (unlike other approaches that consider two nested optimization loops) and the order 
defect problem can be solved without difficulties when the DE method is used. 

4.1 Prescribed timing 

In the problem of path generation n points are given and a mechanism that passes through these points 
is required. In the case of a continuous path a discretization approach can be used. 

When the trajectory is prescribed the input angles for each of the desired points are known, or in the 
worst scenario, are a function of one — the fc-th input angle (typically the first). In such a case there will 
be an extra adjustment parameter, the fc-th or the first input angle. 

The objective function f b for the structural error is: 

n 

/oh = ^||r dl -r gcm || 2 , (4.1) 

^■In fact, it could be possible to reduce the number of parameters by introducing a frame of reference (S'), in which the 
point for the input link joint attached to the fixed link lies on the x' axis, but then all the desired points need to be changed 
to the S' frame; this is not a difficult task. However, it is not necessary to introduce other frameworks. 

2 As far as we know, only [12] addresses this problem. 
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where r di — r goni is the difference between the z-th desired point r di and the i-th generated point r goni 
which can be obtained using Eq. (3.7). In principle, the sum of the length of the geodesies connecting v di 
with r gcn i is what must be minimized, but it is clear that the minimum of Eq. (4.1) implies the minimum 
of the sum of the geodesies connecting the desired and generated points. This can be seen by calculating 
the inner product of the vector = r<& — r gon i with itself, i.e., • ej = 2(1 — cosSi). Where Si is the 
length of the geodesic connecting r di with r gcni . 



4.2 Path generation without prescribed timing 

Unlike the prescribed timing case, in the nonprescribed case the input angles are not known. Therefore, 
it is necessary to introduce as many extra parameters as points to generate. Formally, this does not 
represent a problem, as Eq. (3.7) deals with it, but the optimization turns out to be more complicated, 
because the number of adjustment parameters increases. Let us analyze the situation for the case of n 
desired points. 
Defining 

fob(Oi,02, ■ ■ ■,0n,i3,1,Vi,'P2,<P3,<P4,m,V2,V3,m) = fob(0i,0 2 , ■ ■ .,9 n ,l3,7,(p,ri), 

where (p = (<p 1} ip 2 , <^3, Vi) and rj = (r]i, T] 2 , m)- 
The objective function is then defined as 

n 

fob(0i,0 2 , ■ ■ ■ ,0 n , /3,i, (p,rj) =y~]\\ (6 k ,f3,J,(f,T])\\ 2 . (4.2) 

k=l 

So, the design variable vector is 

X D = {0i, 2 , . . .,0 n ,(3,'y,(p 1 ,(p 2 ,<p3,tp4,rii,r) 2 ,ri 3 ,ri4} ■ (4.3) 

The domain for each of the adjustment parameters in the design variable vector is given as follows 

6 k E [0,2tt); G [0,2^); 7 e [-ir,ir);<p k E [0,2^); V E [0,tt]. 

The next step is to minimize the objective function given by the expression (4.2). Care must be 
taken due to the order defect problem. It would be possible to find a mechanism that passes through the 
desired points, but not necessarily in the required order. Order is imposed because the input link has an 
increasing rotation. This problem is solved by requiring the input angles to be in ascending order as: 

01 < 02 • • • < 0n- (4.4) 

Now, if an evolutionary algorithm is used to solve the optimization problem, all the variables of the 
vector X.d must be randomly chosen. Then it is likely that the 8 angles do not result in ascending 
order. To address this, we could penalize the f b function (which is a common approach in dealing with 
constraints [17]). However, if n is relatively big the penalization approach would not be appropriate, since 
for practical purposes all individuals would be penalized (the probability of finding one that would not, 
is 1/n!). In other words, the evolutionary algorithm would not have individuals to evolve. 

A better approach may be found in [18], which consists of discretizing the search space for the angles, 
and then looking for the first angle in the first partition and so on. However, there is no guarantee that 
the minimum will lie in the discretized space, as it might happen that two of the best fitting angles are 
in the same partition. Since each angle is searched in a different partition the minimum would never be 
found. Another approach — which is used here — is to apply the DE method (as in [19] ) to manipulate 
individuals in such a way that satisfies the constraint (4.4). 
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5 The Differential Evolution Algorithm 



Algorithms that are inspired by natural evolution are known as evolutionary algorithms. In such methods 
there is a population which is susceptible of mutation, crossover and selection. 

The main applications of the evolutionary algorithm have been to optimize functions. The optimizing 
function does not need to be a diffcrcntiable function and the state space of possible solutions can be 
disjoint and can encompass infeasible regions [20]. This is an advantage since in the synthesis of the 
spherical mechanism the discriminant in Eq. (3.4) can become negative. An individual that make the 
discriminant negative can not be a mechanism that optimize the objective function, then the evolutionary 
method will start to look for a better individual in another region of the domain. Among the evolutionary 
algorithms, DE is known for its simplicity and for the excellent results it allows. Below, the original 
version of the method is outlined [17]. 

1. The population is described by: 

Px, g = (xj, s ), i = 1,—m; g = 0, ...g max 
x, : , - (x{. g ), j = l, ...£>; (5.1) 

where D, to and g max represent the dimensionality of x, the number of individuals and the number 
of generations respectively. 

2. Initialization of population: 

xi fi = rand\0,l)-(b^-b{) + b{. 

Vectors by and b L are the parameter limits and randP(0, 1) is a random number in [0, 1) generated 
for each parameter. 

3. Mutation: 

v i;g — x r ;g + F • (x ri;g — X r2;g ). (5-2) 

x rn;9 is called the base vector which is perturbed by the difference of other two vectors. 
r o, r i, r 2 € {1, 2, ...to}, n ^ T2 ^ ^3 7^ i . F is a scale factor greater than zero. 

4. Crossover: 

A dual recombination of vectors is used to generate the trial vector: 

„ _ „,j _ J v ig ii(rand 3 (0, 1) sC Cr or j = j rand ) 

_ U ^ ~ \ 4 9 otherwise. ^ 

The crossover probability, Cr € [0, 1], is a user-defined value, j r and € [1, D]. 

5. Selection: 

The selection is made according to 

_ f U lig if /(U; ;9 ) S$ f(Xi, g ) , 

Xl;9+1 _ \ x,. 3 otherwise. { * > 

The method just described is known as DE/rand/l/bin. There are variants of it. For example, when 
F is chosen to be a random number the variant is called dither. 
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Tabic 1: Desired points for the prescribed timing path generation example. 



Point number Point Point number Point 



1 


(0 


85737, 


-0 


18481, 





48037) 


33 


(0.7887, 


-0.60370, 0.11578) 


2 


(o 


82985, 


-0 


20167, 





52030) 


34 


(0.80152 


-0.59270, 0.07900) 


3 


(o 


80241, 


-0 


21996, 





55478) 


35 


(0.81378 


-0.57959, 0.04311) 


4 


(o 


77567, 


-0 


23967, 





58389) 


36 


(0.82552 


-0.56433, 0.00841) 


5 


(o 


75011, 


-0 


26056, 





60785) 


37 


(0.83678, 


-0.54700, -0.02478) 


6 


(o 


72607, 


-0 


28244, 





62693) 


38 


(0.84759, 


-0.52763, -0.05611) 


7 


(o 


70381, 


-0 


30515, 





64152) 


39 


(0.85807, 


-0.50641, -0.08530) 


8 


(0 


68352, 


-0 


32833, 





65193) 


40 


(0.86819, 


-0.48344, -0.11200) 


9 


(o 


66533, 


-0 


35185, 





65844) 


41 


(0.87804, 


-0.45889, -0.13596) 


10 


(o 


64933, 


-0 


37537, 





66144) 


42 


(0.88763, 


-0.43304, -0.15689) 


11 


(o 


63559, 


-0 


39867, 





66115) 


43 


(0.89704, 


-0.40611, -0.17448) 


12 


(o 


62415, 


-0 


42159, 





65781) 


44 


(0.90626, 


-0.37837, -0.18848) 


13 


(o 


61504, 


-0 


44389, 





65167) 


45 


(0.91537, 


-0.35022, -0.19867) 


14 


(o 


60833, 


-0 


46541, 





64293) 


46 


(0.92433, 


-0.32193, -0.20481) 


15 


(o 


60396, 


-0 


48596, 





63170) 


47 


(0.93322, 


-0.29396, -0.20667) 


16 


(o 


60196, 


-0 


50548, 





61819) 


48 


(0.94196, 


-0.26667, -0.20400) 


17 


(o 


60230, 


-0 


52381, 





60241) 


49 


(0.95052, 


-0.24048, -0.19667) 


18 


(o 


60485, 


-0 


54085, 





58448) 


50 


(0.95885, 


-0.21581, -0.18441) 


19 


(o 


60959, 


-0 


55656, 





56448) 


51 


(0.96685, 


-0.19315, -0.16707) 


20 


(o 


61637, 


-0 


57081, 





54244) 


52 


(0.97430, 


-0.17289, -0.14452) 




/ n 

(U 


6ZO00, 


n 
-U 


O0060, 


n 



K1 QA A \ 

0lo44J 


fro 
OZ 


(^0.98096, 


-0.10041, -U.11606J 


22 


(o 


63530, 


-0 


59489, 





49248) 


54 


(0.98652, 


-0.14104, -0.08315) 


23 


(o 


64700, 


-0 


60456, 





46467) 


55 


(0.99052, 


-0.13007, -0.04430) 


24 


(o 


65989, 


-0 


61259, 





43507) 


56 


(0.99244, 


-0.12263, -0.00019) 


25 


(o 


67370, 


-0 


61896, 





40378) 


57 


(0.99174 


-0.11870, 0.04874) 


26 


(o 


68811, 


-0 


62363, 





37093) 


58 


(0.98774 


-0.11822, 0.10185) 


27 


(o 


70293, 


-0 


62652, 





33674) 


59 


(0.98000 


-0.12085, 0.15807) 


28 


(o 


71789, 


-0 


62759, 





30133) 


60 


(0.96819 


-0.12626, 0.21604) 


29 


(o 


73274, 


-0 


62678, 





26500) 


61 


(0.95226 


-0.13415, 0.27426) 


30 


(o 


74737, 


-0 


62404, 





22800) 


62 


(0.93252 


-0.14411, 0.33115) 


31 


(o 


76167, 


-0 


61933, 





19059) 


63 


(0.90956 


-0.15600, 0.38515) 


32 


(o 


77544, 


-0 


61256, 





15307) 


64 


(0.88422 


-0.16959, 0.43519) 



6 Example 1. Prescribed timing path generation 

The example presented in this section consists in finding the spherical four-bar mechanism that passes by 
sixty four prescribed points. This problem has been studied originally in [11], and later in [13] when the 
desired points lie on a unit sphere. From [13] we take the data to be used, which we reproduce in Table 
1 for quick reference. In the aforementioned references, the problem is solved using the atlas method. 

In order to optimize the objective function we use the DE/rand/l/bin method with the dither variant. 

The first angle Q\ is not known, therefore it will be an adjustment parameter. The other angles are 
given by 

2-k 

0* = 0i + — (fc-l);fc = 2,3,...,64. (6.1) 

The design variable vector is 

X = {9 1 ,f3,j,(pi,ip 2 ,(p3,f4,Vi,V2,V3,Vi} ■ (6.2) 

The objective function is constructed as explained in Section 4.1. The optimal value for the objective 
function that was obtained by the DE method is fob = 3.3 x 10~ 8 . In Table 2, we show the optimal 
design variables. The DE method was coded in Fortran, the population size was of 100 individuals and 
a number of 10 000 generations were used. Fig. 3 shows the desired points, the generated trajectory and 
the obtained mechanism. 
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Table 2: Parameter values of the design variable vector for 64 points with prescribed timing. 



0i 


fl 7 


fi 


V?2 


<P3 


ip 4 


m 




% 


V4 


0.48867 


0.23066 0.47437 


0.000009 


0.38828 


0.19646 


0.97780 


1.57081 


1.46619 


0.66128 


1.34474 



Table 3: Lengths of the links for the example of 64 points with prescribed timing. 



Crank link 


Coupler link 


Oscillator link 


Frame link 


0.40142 


0.82033 


0.92503 


0.99484 



6.1 The Grashof Condition 

Although the Grashof condition [21-23] can be implemented by a penalty term in the DE method, wc 
decided to let the method evolve freely. The Grashof condition is verified at the end of the execution 
program with the purpose of ensuring the mobility of the mechanism in the sense of [22,23]. Using the 
data in Table 2, we can calculate the link lengths. For that purpose, we must find the points for the joints 
of the links (using Eq. (3.8)) and then use the inner product. Table 3 shows the obtained link lengths. 
We can see that the obtained mechanism satisfies the Grashof condition. 

7 Example 2. Path generation without prescribed timing 



Table 4: Angles 9 of the design variable vector for 64 points without prescribed timing. 



0i 


02 


9 3 


04 


05 


06 


07 


08 


2.98039xl0~ 7 


0.0977871 


0.195564 


0.293198 


0.390854 


0.488426 


0.58619 


0.684141 


09 


#10 


flu 


012 


013 


014 


015 


016 


0.782192 


0.880383 


0.978634 


1.07687 


1.17502 


1.27299 


1.37043 


1.4675 


017 


018 


019 


020 


021 


022 


023 


024 


1.5642 


1.66089 


1.7576 


1.85419 


1.95116 


2.04815 


2.14514 


2.2422 


025 


026 


027 


028 


029 


030 


031 


032 


2.33934 


2.4366 


2.53394 


2.63134 


2.72876 


2.82621 


2.92382 


3.02148 


#33 


034 


035 


036 


037 


038 


039 


040 


3.11924 


3.21722 


3.31511 


3.41321 


3.51148 


3.60982 


3.70845 


3.80714 


041 


042 


043 


044 


045 


046 


047 


048 


3.90616 


4.0052 


4.10437 


4.20389 


4.30343 


4.40325 


4.50295 


4.60259 


049 


050 


051 


052 


053 


054 


055 


056 


4.70184 


4.80104 


4.89973 


4.99843 


5.09691 


5.19563 


5.29455 


5.39354 


057 


058 


059 


060 


#61 


062 


063 


064 


5.49265 


5.59183 


5.69105 


5.7901 


5.88913 


5.98799 


6.08645 


6.18493 



As before, we address the problem of finding the spherical four-bar mechanism that passes by sixty-four 
points. But unlike the previous section, we now deal with the non-prescribed timing case where knowledge 
of Eq. (6.1) is not assumed. The purpose of this example is twofold: first, to obtain the mechanism that 
accomplishes the desired trajectory; and second, to see whether the method can realize that the angles 
are equally spaced by the amount of 27r/64. This is a complicated problem that will test the effectiveness 
of the proposed methodology. 

For this case the design variable vector is given by 

X = {01,02, • • • ,064, P, 7, <Pl, <P2,<fi3,V4, m,V2,V3,Vi} ■ (7-1) 
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Figure 3: Desired points, generated trajectory and obtained mechanism. Prescribed path generation. 



Table 5: Remaining parameter values of the design variable vector for 64 points without prescribed timing. 



p 


7 


<Pi 


</>2 


¥>3 


(Pi 


Vi 


V2 


m 




0.293895 


0.455927 


3.10628x10"° 


0.301721 


0.358277 


1.01732 


1.60055 


1.31837 


0.456631 


1.35054 



The construction of the objective function is explained in Section 4.2. As before, the method used 
to solve the optimization problem was DE/rand/l/bin, coded in Fortran with a population of 300 
individuals and a number of 50 000 generations. The order defect problem was handled as in ref [19]. The 
optimal value for the objective function is fob = 5.7 x 10~ 6 . The parameter values for the 9 angles are 
shown in Table 4, the remaining adjustment parameters are shown in table 5. As in the previous example, 
the obtained mechanism satisfies the Grashof condition. 

In order to see if the obtained angles are in agreement with Eq. (6.1) we have calculated the differences 

AO = {9 2 - 9 U 9 3 - 9 2 , . . . , 64 - 9 63 }. (7.2) 

This set of differences have a mean value of /i = 0.0981734. Notice that /i — 27r/64 = 1.3 x 10~ 6 , so the 
error by taking fi instead of 27r/64 is smaller than the accuracy with which the points themselves are 
given. 

We close this section by saying that it is possible to reduce the value of fob by increasing the number 
of generations. However, this is unnecessary since we can realize that the values of Tables 2 and 5 are 
(for practical purposes) very similar 3 . 

8 Conclusions 

By constructing the parametric equations of the spherical geodesies, the generated points of a spherical 
mechanism can be determined in an easy way and the optimum synthesis can be done in a concise manner. 
The article shows in detail how the construction of the objective function for the task of path generation 
without prescribed timing can be done. Using Differential Evolution the structural error was optimized for 

3 Although, we could find another mechanism, since it is well known [24] that the spherical synthesis has not a unique 
solution, in our case we find the same solution. The fact that the design parameters 9\, tpz, and 773 do not have similar 
values for both mechanism, is regardless, since what really matters is that the links lengths for both mechanisms are very 
similar. Also notice that the initial points for the input and output links joints are practically the same for both mechanisms. 
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two examples. The first example considered the prescribed timing path generation for sixty four points. 
The method was capable of finding an optimal value for the objective function of 3.3 x 1CT 8 . In the 
second example, we solved the problem of path generation without prescribed timing for the same sixty 
four points. The non-prescribed timing case demands the introduction of sixty four additional adjustment 
parameters for the angles. Hence a total of seventy four parameters are required to construct the objective 
function. The optimization process was done by minimizing a single objective function without the need 
for separate optimization stages as has been done in previous studies. The method was capable to obtain 
an optimal value of 5.7 x 1CT 6 for the objective function and also to obtain the correct prescription with 
a difference of 10~ 6 from the exact value. 
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